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Abstract 

Using random variables as motivation, this paper presents an exposi- 
tion of formalisms developed in [RT1, RT2] for the classical umbral calcu- 
lus. A variety of examples are presented, culminating in several descrip- 
tions of sequences of binomial type in terms of umbral polynomials. 
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1 Introduction 

The system of calculation now known as the "umbral calculus" originated with 
Blissard in the nineteenth century in informal calculations involving the "low- 
ering" and "raising" of exponents. The work of Rota and his collaborators 



in [MR, RKO, RmR] and other works formalized these methods in the modern 
language of linear operators and Hopf algebras. While this clarified the under- 
lying theory, it rendered the original nineteenth century work no more easy to 
read or check. In |RT1| , |RT2| , the original classical notation was revived and 
extended so as to be rigorous by modern standards. In [RST|, the first attempts 
were made to apply this newly revived classical umbral calculus to one of the 
most significant successes of the modern theory, namely the study of sequences 
of binomial type. 

The purposes of this paper are twofold. Since much of this paper is expos- 
itory, no prior knowledge of umbral calculus in any of its guises is assumed. 
To start, we develop the modern formulation of the classical umbral calculus 
in analogy w ith the idea of a random variable. This renders the definitions 
of RT1 , RT2] transparent. In the second part of this paper, we introduce a new 
operation on umbrae arising naturally from the analogy to random variables. 
We show that all sequences of binomial type and all umbral maps arise directly 
from the application of this operation. We further apply the tools of classical 
umbral calculus developed in the first part of the paper to provide several other 
compact presentations for sequences of binomial type. 

The author has attempted to document at least the recent history of the 
main results and definitions contained herein. 
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2 Random variables and the classical umbral cal- 
culus 



Fundamental to the classical umbral calculus is the idea of associating a se- 
quences of numbers do, a±, a%, . . . to an "umbral variable" a which is said to 
represent the sequence. To be slightly more formal, the umbral calculus re- 
lies on associating the sequence ao, <zi, ai, • • • to the sequence 1, a, a 2 , a 3 , . . . of 
powers of a. 

This kind of association is familiar in modern mathematics: To any random 
variable G, we associate a sequence of numbers 1, gi, . . . where gi is the ?th 
moment of G. Specifically, we are defining the sequence 1, (71,32, ■ ■ ■ by apply- 
ing the expectation operator, E, componentwise to the sequence 1,G, G 2 , . . . 
consisting of powers of the random variable G. 

We will proceed to carry this analogy further, calculating with random 
variables in precisely the way we will later be using umbral variables. We 
start by letting G be a random variable distributed uniformly over the interval 
[0, 1]. The sequence of moments 1,91, g%, ■ ■ ■ associated to G is thus given by 



If we let p(t) £ C[t] be any polynomial with complex coefficients, it is imme- 
diate that 



where A is the forward difference operator Ap(t) = p(t + 1) — pit). A comment 
on notation: since A is defined as an operator on the ring of polynomials in t, 
Ap(0) can only be interpreted as (^pj (0) or as 0. We adopt the former reading. 

Since the calculation in equation [j] only required that p(t) was differentiable, 
it could just as well have been carried out for a polynomial with coefficients in 
some other integral domain or indeed for a polynomial in t whose coefficients 
contained various random variables which were independent of G. So suppose G' 
is a random variable independent of and identically distributed to G. Consider 
p(G+G r ) as a polynomial q(G) in G with coefficients in C [£?']. Using equation |l], 
the expected value, averaging over values of G, of p"(G+G') — q"{G) is Aq(0) = 
Ap'(G'). Applying equation [j] again, recalling that the derivative D and A = 
e D — I commute and that G' is identically distributed to G, gives E[Ap'(G')] = 
A 2 p(0). Somewhat more suggestively, this calculation can be written 



The first property to observe here is that the independence of G and G' 
matters, as 



9n = Jo9 n d9= ^tt- 




(1) 



E[p"(G + G')] = E[Ap'(G')] = E[Ap(0)] . 



(2) 



E[p"(G + G)] = E[Ap"(2G)] = 




since p"(2t) 



Dt p 2 . where D t is the derivative with respect to t. 
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The second important property is the rather trivial observation that we 
can calculate the expectation of a polynomial in several random variables all 
independent and identically distributed to G simply by knowing the moments 
of G. For example, by equation ^, we know that E[20(G + G') 3 ] = A 2 i 5 | t=0 . 
This could be evaluated directly. Alternatively, and denoting E[G k ] by G&, it 
could be evaluated as 

20 • E[G 3 + 3G 2 G' + 3GG' 2 + G' 3 ] = 20(G 3 + 3G 2 Gi + 3GiG 2 + G 3 ) 

= 20(2G 3 + 6G 2 Gi) 




Applying the above observations quickly reconstructs the moment generating 
function E[e Gz ] for G. Since D t sJl = e tz , equation |l| implies E[e Gz ] — e '~ e . 

For the duration of the next calculation, we are going to make some assump- 
tions that simply do not hold within the confines of probability theory. The 
remainder of this section will be devoted to describing how to replace random 
variables with "umbral variables" in a way that makes the following calcula- 
tions legitimate. So, for the moment, let us assume that there is an object B 
that behaves much like a random variable. Let us call this object an "umbral 
variable." We treat it just like a random variable, but stipulate both that it is 
independent of B and that B + G = 0. With these stipulations, we find that by 
independence 

E j e (G+s) 2 = ft[e Gz e Bz ] = E[e Gz ] E[e B2 ] = - E[e Bz ] . 

But since G + B = 0, the left-hand side above is just 1. We have just cal- 
culated that E[e Bz ] = -tzT' Since this is the exponential generating function, 

Sfe>o ^ fc Tp f° r ^ ne Bernoulli numbers, we find that if the above calculation 
can be made rigorous, then E[_B fc ] = B^. 



Since calculations such as the previous are too useful to abandon (see [RT2 
for a variety of examples involving the Bernoulli numbers) we define umbral 
variables or umbrae which formalize the roles of both the random variable G 
and the new object B in the preceding calculation. 

Just as a random variable is usually capitalized, we will typically distinguish 
our umbrae by writing them as Greek letters, e.g. a, (3, a'" , .... Let us denote the 



collection of whichever umbral variables we will be using by A. See [RT1, RT2 
for the relevant, and straightforward, technical details. In practice, when we 
introduce a new umbra, say a, we specify explicitly or implicitly how E acts on 
it, namely what values E[a fe ] takes for each k. Any two distinct umbrae in A, 
say a and 7 or a and a' will acts as independent random variables, regardless of 
how E acts on them. Generalizing, any collection of distinct umbrae will behave 
as independent random variables. 

Formally, this can be accomplished by defining a linear evaluation map, 
E : F[A] — > F on the polynomials F[A] in the umbrae with coefficients in a 
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suitably chosen commutative ring F. We require that E is F-linear, E[l] = 1, 
and E[M • M'] = E[M] E[M'] for any two monomials, M and M' in CA such 
that no umbra appears to nonzero power in both M and M'. This map was 
called eval in |RT1| , |RT2| . 



We call p, q € F[„4] umbrally equivalent, written p ~ q when E[p] = E[g]. 
Analogous to the notion of identically distributed random variables, we define 
p, q G F[„4] to be umbrally exchangeable, when p k ~ q k for all k > 0. If for 
example we have a = 3, then E[a fe ] = 3 fe ; this is consistent, in the analogy 
to random variables, with considering a analogous to a random variable which 
always takes on the value 3. We not that equality implies exchangeability which 
implies umbral equivalence. The converses are false. 

We define p,q G F[_4] to be independent when no umbra appears in both 
p and q. More formally, an umbra that appears to a nonzero power in some 
monomial with nonzero coefficient in p does not appear to a nonzero power in any 
monomial with nonzero coefficient in q. For example, a 2 +aa' and /3(3' 2 — /3+a'" 
are independent, but a 2 + aa' and /3/3' 2 — /3+a are not independent. Nor are the 
falling factorials (a)( n ) and (a)(„_i) independent. On the other hand, at 2 + j3 
and 7 — t are independent where t is some element of F. 

In random variables, we can usually substitute one identically distributed 
random variable for another (modulo independence constraints). Similarly, we 
can substitute exchangeable umbrae as per the following lemma. 

Lemma 1 (RT1) If a polynomial p(t) is independent of two exchangeable um- 
brae a and a', then p(a) =p(a'). D 

This substitution lemma holds equally well if a or a' is replaced by an umbral 
polynomial p G F[.4]. 

For the duration of this paper, we let e be an umbra such that e k c± Sq^ 
where S in the Kronecker delta. As long as we work with polynomials in the 
umbrae, there is no harm in defining =1. Under this convention we consider 
= e. This is consistent with the convention that 1 is the 0th moment of a 
random variable which always takes on value 0. 

To pick up our earlier example, let 7 be an umbra such that 7* ~ ^xj and 
let (3 be an umbra such that 7 + [3 = 0. It is an easy exercise to see that given 
any umbra 7 such an inverse umbra can be found recursively. Here we have 
made formal in 7 and (3 exactly the properties we had assumed for G and B. 
We have 7 and (3 are independent and (3 + 7 = 0. 

Extending our notions of independence, exchangeability, and equivalence 
coefficientwise to formal power series (see Q for a general but technical treat- 
ment) we can duplicate the computation we did for G and B. By the substitution 
lemma we have e^ +7 ^ ~ e°' z = 1. Thus we have e' 3z e 7Z ~ 1. By independence, 
this implies that e^ z ^—^- ~ 1 and hence (by linearity) that e' 32 ~ —rzj- From 
a technical viewpoint, there are a number of ways to justify the first step in the 
preceding sentence. The most direct solution is to apply the substitution lemma 
coefficientwise to the formal power series in z. A general approach which views 
multiplication by e 7Z z z _, as a linear operator "equivalent" to the identity is 
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given in PJ. The intuition behind both these proofs is that since e 7Z and e@ z 
are independent, E can be applied in two stages first to 7 and then to /3, anal- 
ogously to finding the expectation by first averaging over one random variable 
and then over another independent random variable. 

As a demonstration of these techniques we rewrite in modern umbral nota- 
tion the first example in jBlj , one in the series of papers in which Blissard during 
the 1860's introduced his "representative notation" — the umbral calculus. To 
point out just how closely the modern language captures Blissard's 19th century 
original, we present most of this example in Blissard's own words. 

Blissard starts with the problem "Required to expand { i og (i+ x ) } 

He then lets "| log{ x 1+x) } = 1 + P\x + P 2 x 2 + ■■■ + P n x n + &c." and 

defines U n to be the coefficient of ^- in I e J where is an ordinary 
variable. He observes that 

e 6 — 1 \ m 9 2 6 n 

— = i + u 1 e + u 2 — + --- + u n — + &c. 

a J 1-2 1 ■ 2 ■ ■ ■ n 

= e ue (by representative notation) . " 

In modern language, he is letting U be an umbra such that U n ~ U n , 
also his "=" would be replaced with The next operation takes 

place purely on the level of formal power series. Blissard substitutes 
log(l -fx) for 6 and finds that 

n m 

= (l + x) u " 



log(l 



where again the only change necessary to modernize his work is to 
replace "=" with If we "equate coefficients of x n , then P n = 

U(U-l)(U-2y--(U-n+l) „. again we wQuld replace U= » with 

The preceding formula for P n has the advantage of being extremely compact. 
Blissard concludes with an expansion of it, and we shall proceed likewise, though 
our precise techniques are somewhat more umbral than those Blissard used. 

With 7 as before, we have ^ ~ e ~< e . Thus (^"j™ ^ e (i' +n" +~W")0 

where 7' + 7" + • • • + 7"' is a sum of m distinct (and thus independent) umbrae 
each exchangeable with 7. We conclude that U = 7' + 7" + • • • + 7"' and thus 
P n ~ ( 7 +7 + n " +1 ) • We conclude with the following formula for evaluating the 
powers U n . Since jj r t n t m+n = (m-f n)( m )t", we have, generalizing calculation ||, 
that 



m rim+n 



(m + n) [m) U n ~ (m + n) (m) ( 7 ' + 7" + • • • + 7"')" * A™ 

This last is better known as ml S(m + n, m), where S(n, k) is the Stirling number 
of the second kind counting the number of set partitions of an n-set into k parts. 
So U n ~ S(m + n, m) / ( m r ^™) • We can use this, together with the expansion of 
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the falling factorials in terms of Stirling numbers of the first kind to derive 

'•' + ' :„ + y ") = h t »(»■ »)(V + - + vr = 5 1 *"S U| ■ 

' ' fc=0 ' fc=0 \ m / 

3 Umbral presentations of Appell sequences 

Historically, the objects of interest in umbral computations were of course se- 
quences of numbers or polynomials. For our present purposes, this means that 
we will primarily be studying the "moments" -E[a fc ] of an umbra a. We say that 
the umbra a represents a sequence do, cti, a^, ■ ■ ., a; £ F, when a k ~ au for all 
integers k > 0. Necessarily this implies that ao = 1. An umbral presentation of 
a sequence oi, &2, . . . of elements in F is any sequence qi, q%, . . . of polynomials 
in F[A] such that qi ~ dj for i > 0. Throughout this paper we freely assume 
that, for any sequence in F, we can find infinitely many umbrae representing 
the given sequence. 

Now let F be \c[x, y] where k is a commutative ring containing Q. The 
remainder of this paper will focus on umbral presentations for sequences of 
polynomials. For example for any umbra a, we can define a sequence of poly- 
nomials s n (x), n — 0, 1,2, . . ., by s n (x) ~ (x + a) n . This definition immediately 
yields the calculation 



s n (y + x)c(y + x + a) n = ]T (^jy\x + a) n - 1 ~ Yl (") 



y l s n -i{x). (3) 



A sequence of polynomials so(x), si(x), S2[x), . . . is said to be an Appell 
sequence when it satisfies the identity, 



n (y + x)=J2( n )y is n-i(x) (4) 



given by equation || for all n > 0. We shall call an Appell sequence s n (x) 
normalized when S\(x) is monic. Any Appell sequence may be rewritten as a 
normalized Appell sequence by replacing s n (x) with s n (x)/s' 1 (0). Here, s[(x) is 
the first derivative of s(x). We hold with this notation for derivatives throughout 
this paper. In the literature, Appell sequences are frequently defined to be 
normalized. 



Proposition 2 ( flFtT2j ]) A sequence, Sq(x), Si(x), . . ., of polynomials in k[x] 



with s n (x) having degree n is a normalized Appell sequence iff there exists an 
umbra a such that s n ~ (x + a) n for n > 0. 

Proof. The if direction is given by calculation ^|. 

(only if): Replacing x with in the defining equation ^ shows that in an 
Appell sequence, each polynomial s n (y) can be recovered from the sequence 
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of values so(0), Si(0), . . .. Choosing an umbra a that represents this sequence 
guarantees (x + a) n ~ s n (x). □ 

Similarly we have the standard result that a sequence sq(x), si(x), . . . of 
polynomials, s n (x) of degree n, is an Appell sequence iff 

s' n (x) = n ■ 8 n -!(x) (5) 

for all n > 0. The only if direction follows since D x (x + a) n = n(x + a)™ -1 . To 
show the if direction, we observe that any sequence of polynomials satisfying 
equation^ is determined by the sequence of values of sq(0), si(0), sa(0)> ■ ■ ■ an d 
apply the argument in the preceding proof. 

The sequences of polynomials with which we will most concerned in this 
paper are those of "binomial type," i.e. sequences of polynomials which satisfy 
an analog of the binomial theorem. Before approaching this topic however, we 
need to lift another tool from random variables to umbral calculus. 



4 Sums of random variables and the "dot" op- 
eration on umbrae 

Suppose that X is a random variable. If n is a positive integer, one can of 
course run n trials of X and sum the results. Denote the sum by a new random 
variable n.X. Thus n.X has the same distribution as X\ +X2 + ■ ■ - + X n where 



the Xi are all independent and identically distributed to X. In [RT1, RT2| the 
corresponding notion n.a was defined for an arbitrary umbra a. In particular njx 

is itself an umbra and it is defined to be exchangeable with ax + a.% H V a n 

where a, = a for each i. Similarly, for any umbral polynomial p G F[A] we 
define a new umbra n.p which is exchangeable with 71.7 where 7 is any umbra 
satisfying 7 = p. It is worth emphasizing that n.p is itself an umbra. Thus a 
5. a, 3. a, and 5. (a + /3) are all distinct (hence all independent.) It is however 
clear from the definitions that 5. (a + 0) = 5. a + 5.(3. 



We recall a technical consideration from [ RT1 , RT2 1 . The set of all umbrae 
A will be decomposed as a disjoint union A = Ao W A\ , umbrae in .4.1 are called 
auxiliary umbrae. Umbrae of the form n.p are auxiliary umbrae. This detail 
will be given more attention below. 

It is an easy observation that if g(z) ~ e az (this object is analogous to the 
moment generating function of a random variable) then 

e n. a z = e ai+-+a„ » = JJ e <*« * ~ g ( z )» ; 

i=l 

where the last equivalence uses the independence of the e° iZ, s. 
Even more directly, we see that 

(m + n).a = m.a + n.a' (6) 
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where a! = a. As a consequence, if for each positive integer n we define a 
sequence fo(n), fi{n), f^in), ■ . ., by fi(n) ~ (n.a) 1 , then equation O implies that 



((m + n).a) k ~ (m.a + n.a') k = ^ y. j (m.a)\n.a') 



i\k—i 



hence 



A(m + n)= 2 (7) /*-<(»)■ (7) 



This kind of generalized binomial theorem will be explored further in the next 
section. 

By way of introduction to the first new definition of this paper, we consider 
the following generalization of nji . Let X be some random variable and let Y be 
a random variable which only takes positive integer values. Run one trial of Y, 
then run Y trials of X and sum the results. We define Y.X to be a new random 
variable whose distribution is identical to Xi+X 2 + - ■ -+Xy', for convenience, we 
are defining the X^s to be independent random variables identically distributed 
to X. Observe that if Z is another random variable taking only takes positive 
integer values, then according the preceding definition (Y + Z).X has the same 
distribution as Xx+X 2 + - ■ -+Xy+Xy+i + - ■ -+Xy+z an d hence as Y^i+ZX 2 - 
On the other hand, X ,(Y + Z) does not in general have the same distribution 
as Xi.Y + X 2 .Z\ in the first expression only one trial of X is made and in the 
second one makes two trials of X. Of course if X always returns the same value, 
say X = n, this causes no trouble and n,(X + Y) is identically distributed to 
n.X + n.Y. 

Just as the definition n.X for random variables generalizes to n.j for um- 
brae, we would like a generalization of the random variable X ,Y to umbrae. 
This generalization should satisfy results analogous to those observed above for 
random variables. 

The generalization relies on a simple observation first applied to the umbral 



calculus by Nigel Ray in [Ra2] 



Proposition 3 J/7 G Aq is an umbra, and n is a positive integer, then 
E[(n.~f) k ] is a polynomial in n. 

Proof. This is equivalent to the observation that if g(z) is in k[[z]], the ring of 
formal power series in z, and if g(z) ~ e 72 , then e n ' lz ~ g(z) n — e nios ^ z '' and 
the coefficient of z k /k\ in the last is a polynomial in n. □ 



Alternately, we could have observed that if 7* ~ dj, and each 7, = 7 for 
i = 1 . . . n then a^a 1 ^ ■ ■ ■ a 1 ^ appears in E[(7i + • • • + 7n) fc ] as many times as 
there are monomials in the expansion of (71 + ■ ■ ■ + "f n ) k containing exactly ij 
jth powers. But this says that (n.~f) k is umbrally equivalent to 

i\ times iu times 
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where i = n- (ii + . . . + i k ) and („_( il+ ...+ ijfc ) jil) ... )ifc ) is a polynomial of degree 
(ii + . . . + i k ) in n. 

But this is precisely what we need to make sense of replacing n with a. 

Definition 4 Let ajE Aq be umbrae, and define g(z) G F[[z]] by g(z) ~ e JZ . 
Let q lt k{n) be the coefficient of z k /k\ in e ni °^ 9<yZ ^ . Define a.j G A\ to be a 
new auxiliary umbra, such that (a.j) k ~ q^^(a). 

In general, ifp,q G F[.4o] a^e umbral polynomials, we define an auxiliary umbra 
p.q G Ai by p.q = a. (3 where a = p and f3 = q. 

Equivalently we could have defined (a. 7) by replacing n with a in equation |8[ 
The definition immediately implies that e a ' lz ~ e Qlog ( ff ( z )' = (g(z)) a . Sim- 
ilarly, if a(z) G k[[z]] is defined by a(z) ~ e az , then e a ' lz ~ a( log(<7(z))) . It 
is a straightforward exercise in probability theory to show that if a(z) is the 
moment generating function of a random variable A taking only positive integer 
values and if g(z) is the moment generating function of a random variable B 
then B\ + • • • + Ba also has moment generating function a( log(<?(z))) . 

It follows that under this definition O.7 = e = which is what one would 
expect from the analogy to random variables. 

We now state the promised analogues to the standard results on random 
variables. 

Proposition 5 Letp,q,r G F[_4] be umbral polynomials. Ifp,q are independent 
then (p + q).r = p.r + q.r. 

Proof. By definition, and the substitution lemma (Lemma it suffices to prove 
that (a + /3).7 = a. 7 + (3,-f for any distinct umbrae a,/?, 7, i.e. that the fcth 
powers of each side of the displayed equation are umbrally equivalent for all 
k > 0. Letting q ly k{n) be the polynomials from Definition^, it suffices to show, 
for all k > 0, that the identity q lt k{ot + P) — Yli {j) c li,ii. a ) c li,k-i{P) holds purely 
on the level of polynomials in variables a,j. But this follows since equation]^ 
says this identity holds with a, 7 replaced by any pair of positive integers. D 



As remarked above, we cannot expect that p. (q + r) = p.q + p.r will hold 
in general. However, we record the special case where p involves no umbrae. 

Proposition 6 Let a be an element o/F. Let q, r G F[.4i] be umbral polyno- 
mials. If q,r are independent, then a.(q + r) = a.q + a.r. 

Proof. The result holds when a is any integer. Repeating the argument in the 
proof of Proposition || shows the identity holds when interpreted in terms of 
polynomials in a. □ 

The importance of independence is illuminated if we examine what fails on 
replacing a with a. and trying to prove that a. ((3 +7) = a. f3 + a. 7. Staying 
with the notation introduced in Definition ^, we would need to show that 
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a. This is not an equality. It only worked when a = a for 
a € F because the substitution lemma tells that qp+ lt k{o>) — 9/3+7, fc (a) and 

Si (i)(?/8,»(a)g7,fc-»(a0 - Si (i)^,»(°)37,fc-*(o) which is equal to g^+^a). 

The special case of a. 7 where a = —n where n is a positive integer is Ray's 
definition in [Fta2| of "negative umbral integers." Let n be a positive integer. 
Since — n.7 + 71.7 = (— n + n).j = e, the umbra — n.7 defined as above is 
precisely the same as the umbral — n.7 defined in RT2fl . 

The same techniques used in the preceding propositions proves the following. 



Proposition 7 Let a,c be in F. If p 6 F[y4] t/ien a.(cp) = c(a.p). 

This points out that —l.a is noi in general exchangeable with —a. The latter 
is exchangeable with —1(1. a). 

Our definition oip.q does not allow for p or q to contain auxiliary umbrae. 
Nevertheless, we would like to be able to manipulate expressions that resemble 
a,(/3.7). Before we extend the notion of an auxiliary umbra to handle this kind 
of construct, we prove the following associativity result. 

Proposition 8 Let a, (3, 7 be umbrae. Define an umbra p by p = a. (3 and an 
umbra a by a = f3."f. We have p."/ = a. a. 

Before presenting the proof, which is a quick calculation with generating func- 
tions, we interpret the result probabilistically. A.(B.C) can be viewed as finding 
A, then running A trials of B.C, i.e. A times we run a trial of B and following 
each trial of B we run that many trials of C . Then we add up all the trials of 
C. In this interpretation (A.B).C differs only in that we run A trials of B and 
then run all the trials of C at once. We could extend this umbrae a, (3, 7 by 
viewing each side as a polynomial in the variables E[a J ], E[/3 l ], E[7*]. 

Alternately, we argue as follows. 
Proof. It suffices to check that e p * 7 z ~ e a ' a ' z . If a(z) , b(z) , c(z) G k[[z]] are given 
by a(z) ~ e az , b(z) ~ e^ 2 , and c(z) ~ e 72 , then this amounts to observing that 
each side is umbrally equivalent to the composition a(z) olog(z) o&olog(z) oc(z). 

□ 

With this lemma in hand, the following definition makes sense. 

Definition 9 Given umbral polynomials P\,...,p n £ k(y4o], inductively define 
the auxiliary umbra pi.p 2 .. ■ -.p n 6 A\ by p lm p 2 .- ■ -.p n = Pup where p = p 2 .- ■ -.p n - 

5 Presentations for sequences of binomial type 

5.1 Sequences of binomial type and sums of umbrae 

The notion of a sequence of binomial type is a direct generalization of equation 0. 
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Definition 10 A sequence of polynomials po(x),p±(x),p2{x), .. . with p n (x) of 
degree n is of binomial type when it satisfies 

Pk(x + y) = ^2\ .)Pi{x)p k -i(y) . (9) 



Such a sequence is normalized when pi(x) is monic (equivalently p\{x) — x). 

Equation ^ arose directly as the umbral expansion of the identity (equation ^ 
that (m+n)ja = mjoe+rw>' . Recall that, for any umbra 7 and any element x G F, 
E[(x.j) n ] is a polynomial in x and that (x + y)«7 = x.j + y.7 where y is also in 
F. As in the proof of proposition |^, raising both sides of the preceding equality 
to the nth power and applying E implies the if direction of the following. 

Theorem 11 Let po(x),pi(x),p2(x) ) . . . be a sequence of polynomials where 
p n (x) has degree n. This is as sequence of binomial type iff it is umbrally 
represented by x."f for some umbra 7. 

Proof. By the remarks preceding the theorem, it suffices to show that any 
sequence of binomial type can be so represented. By standard results, which are 
briefly sketched below, it suffices to show that choosing 7 appropriately allows us 
to choose the sequence D x E[(x.'j) 1 }\ x=0 , D x E[(xsf) 2 ] | x=0 , D x E[(x.j) 3 ] \ x=0 
arbitrarily. 

It is enough to observe that equation ^ tells us that the coefficient of x in 
(x.'y) k is 7 fc + R where R is depends only on k and E[7], . . . , E[7 fc_1 ]. □ 

For completeness, we sketch the fact that knowing p[(0) , p' 2 (0) , ■ ■ ■ deter- 
mines a sequence po(x) , pi(x) , P2{x) , ... of binomial type. Replacing y with 
in equation ^| and recalling that, by degree considerations the Pi(x) are linearly 
independent, which tells us that po{0) = 1 and that Pi(0) = for i > 0. Tak- 
ing the derivative of equation ^ with respect to y and and setting y to gives 
Pk( x ) = J2iZn Q)Pii x )Pk-i( Q )- Since = <*fc,o tnis determines p k (x). 



Following [ RKOfl , the umbral composition a(h(x)), of two polynomial se- 



quences ao(x), ai(x), . . . and b (x), bi{x), ... is the the sequence T(a(0)), T(a(l)), . 
where T : k[x] —> k[x] is the linear operator defined by T(x l ) — bi(x) for all i. 
An umbral operator is defined to be a linear operator U : k[x] — * k[x] such that 
the sequence U(l), U (x), U (x 2 ), ... is of binomial type. The following corollaries 
are immediate. 

Corollary 12 A linear operator U : k[x] — > k[x] is an umbral operator iff there 
exists an umbra 7 such that U(r(x)) = E[(r(x)) .7] for all r(x) £ Q[x}. □ 

Corollary 13 Let po,pi, ■ . ■ and qo,qi, ■ ■ ■ be sequences of binomial type repre- 
sented by x.a and x./3 respectively. The umbral composition of these sequences, 
Po(ct),Pi(cj)7 . . . is represented by x.[3.a. □ 



This makes obvious the fact from [FLKO] that the umbral composition of two 



sequences of binomial type is also of binomial type. 
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5.2 Generalized Abel polynomials 

One of the best known sequences of binomial type has as its degree n polynomial 
the Abel polynomial x(x + na)"" 1 where a is a constant. Generalizing a to be 
an arbitrary umbra a and replacing na with n.a yields the following. 



Theorem 14 ([ RSTj ]) Let p n (x) 6 k[x] be a sequence of polynomials with 



Pi(x) = x and p n (x) of degree n. 

The sequence p n (x) is of binomial type iff there exists an umbra a such that 

p n (x) ~ x(x + n.a)™ -1 . 



The proof in [ R.ST | closely parallels the proof that the original Abel polynomials 
are of binomial type. Here we provide a combinatorial proof. 
Proof. To start with, assume that E[a l ] is always an integer and that ai, . . . , a n 
are distinct umbrae all exchangeable with a. We start by interpreting (x + 
oti + • • • + a n ) n ~ l as a generating function for sequences of length n — 1 on 
n + 1 symbols. By the Priifer correspondence (see for example [||) this is a 
generating function for the number of labeled free trees onn+1 vertices where 
each tree is counted with weight x d ° J}; a f l where vertex I has degree di + 1. So 
x(x + «!+■■■ + a„)™ _1 is the generating function for labeled rooted trees on 
n + 1 vertices where the same weight indicates that vertex I has outdegree di . 
This says that E[x(x + a.\ + • • • + an)™" 1 ] is the generating function where the 
coefficient of x k counts the number of labeled trees on n + 1 vertices where the 
root has degree k and each non-root vertex with outdegree i can be colored in 
any of ai ways. Equivalently, E[x(a; + ai + ■ ■ ■ + a„)" _1 ] counts the number of 
planted forests on n vertices where each vertex with outdegree i can be colored 
in any of a.; ways and where each tree in the forest can itself be colored in any 
of x ways. Let's call this structure a (x, a) degree-colored forest on n vertices. 

So counting the number of ways to form a (x + y, a) degree-colored forest on 
n vertex by the number of vertices, i, in the trees which were colored in one of 
the first x ways gives 

(x + y){x + y + n.a) 11 - 1 ~ ^ (")(*)(* + i.a) 4 " 1 • (y)(y + n - i.a)" -4-1 . 

i ^ ' 

This fact for all positive integers x, y, a\, a^, . . . implies equation |^ as a polyno- 
mial identity. 

To see that indeed all normalized sequences of binomial type arise in this fash- 
ion, it suffices, by the remarks after Theorem [H], to observe that the sequence 
P2(0),P3(0), . . . can be chosen arbitrarily. Indeed p' n (0) — (n.a)™" 1 ~ na n ~ 1 +R 
where R is a linear combination of a±, . . . , a n ~\. D 

The interpretation of x(x + n.a) n ~ l as a generating function for colored forests 
was suggested to the author by Nigel Ray. It generalizes the notion of reluctant 



functions developed by Muffin and Rota in [MR] and is closely related to the 



chromatic polynomials in [RaSW|. 

The calculations used in [RST| to prove Theorem [l4| show the following. 
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Proposition 15 For any umbra a and any n > 1 we find that 

D x e- UaD °> (x(x + n.a)^ 1 ) ~ nx(x + (n - l).a) n - 2 . □ 



Corollary 16 A sequence po{x),p\(x), . . . of polynomials, p n (x) of degree n is 
a sequence of binomial type iff there exists a formal power series g(t) € k[[tj] 
with g(Q) = and g'(0) ^ such that g(D x )(p n (x)) = np n _i(x) for all n > 1. 

Proof. If is monic or ,g'(0) = 1, the result follows immediately from the 

preceding proposition and the remarks after Theorem [H]. 

If the sequence is of binomial type and Px(x) = ax, then so is the sequence 
whose nth term is p n (x)/a n ; if g(t) is the series associated to this new sequence, 
then g(t)/a works for the original. The converse follows similarly. □ 

Formal power series of the sort described above are called delta operators and 
the correspondence between them and their associated sequences of binomial 



type was established in [RKO] 



The transfer formula also arises as an immediate corollary. Since we now 
know that the operator Q ~ De~ 1,aDx is associated to the sequence presented 
asp„(x) ~x(x+oti-\ Vct n ), it follows that p n (x) ~ a;e («i+---+ a! «)- D * = 

«l-D, . . . „a n D x ( n — l\ ^ 



xe 

With the definition of a delta operator in hand, we recall the first expan- 



sion theorem from |RKO|. If Q is a delta operator associated to a sequence 
Po(x) ,pi(x), ... of binomial type and if T = f(D x ) for some formal power series 
f(t) e k[[t]], then 

T = J2T Pk (0)T V - (10) 

fc>0 

The usual proof employs the binomial expansion for p n (x + y) to verify the 
identity f{D y ) = Y,k>o f( D v)Pk(y)kT on Pn{x + y) for all n > 0. Setting 
y = gives the desired result. The fact that f(D y ) and the shift operator e xDy 
commute is used freely. 

We have recalled the expansion theorem in order to derive the identity 

0. = £ri(°)TT (11) 

fe>0 

where po(x),pi(x), ... is the sequence of binomial type associated to a delta 
operator, Q. Thus if Q = f(D), we recover the fact that p' k (0) is the coefficient 
of t k jk\ in the power series inverse, under composition, to fit). This 

relationship together with the transfer formula is used to prove the Lagrange 



inversion formula. See [ MR , RST ] for such derivations 
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5.3 Generalized rising factorials 

Our next presentation generalizes the binomial type sequence of rising factorials, 
x(x + 1) . . . (a; + n — 1). More generally, it is well known that the sequence 
p n (x) = x(x + a) . . . {x + (n — l)a) is of binomial type for all constants and that 
its associated delta operator for a ^ is (I — e~ aD )/a. For the rising factorials, 
this is the backwards difference operator f(x) f(x) — f(x — 1). Our proofs 
closely follow those for the usual rising factorials. 

Theorem 17 Let p n {x) E k[x] be a sequence of polynomials with pi(x) = x 
and p n (x) of degree n. 

The sequence p n (x) is of binomial type iff there exists an umbra p, such that 

p„(x) ~ x(x + (ii)(x + m + (i 2 )(x + m + (i 2 + (1,3) ■ ■ ■ (x + Mi + • • • + (ln-l) 

where (ii, . . . , (i n -i are distinct umbrae exchangeable with (i. 

Proof. We start by showing that all such presentations are of binomial type. 
By induction we have that 

(x + y)(x + y + (ii)(x + y + (ii + (i 2 ) ■ ■ ■ {x + y + p,i + . . . + (i n -i) ^ 

y(y + x + (i)(y + x + (i + ^(y + x + (i + (i x + (i 2 ) ■ ■ ■ (y + x + (x + (i\ + . . . + (i n - 2 ) + 
+x(x + y + (jl)(x + y + (jl + fJ,i){x + y + (i + (ii + (i 2 ) ■ ■ ■ (x + y + (i + \i x + . . . + (i n - 2 ) 
^ yPn-i(x+ (y + M)) + xp n _ 1 ((x + (i) +yj 

= i ^Pn-i-iix^iy + (i) + (" i ^Pn-i-tix + (i)pt(y) 

- E [ ■ )Pn-i-i{ x )y{y + v)(y + a* + a*i) ■ ■ • (y + a* + a*i + • • • + m»-i) + 

i ^ ' 

+ E ( H ■ j X ^ X + A t )( x + A* + Ml) ' ' ' + M + A*l + • • • + A*n-2-*)p»(y) 

- E ( 71 1 + ( H ■ )Pn-i( x )Pi(y) 



- it2 (j _ l^jpn-i( x )Pi(y) + Yl (™ i ^Pn-ii^Piiy) 

- E (^jpn-i( x )Pi(y) ■ 



As before, all normalized sequences of binomial type arise in this fashion, 
since the sequence p' 2 (0),P3(0), • • • can be chosen arbitrarily. Observe that 
p' n (0) ~ A*i(a*i + M2) • • ■ (a*i + • • ■ + A*rt-i) — Mi _1 + R where each term of 
R has degree less than n — 1. □ 
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Proposition 18 Let p be an umbra and let jj,2, . . . be distinct umbrae ex- 
changeable with p. Consider the sequence of binomial type p n (x) ~ x(x + 
jti) • • • (x + fi\ + ■ ■ ■ + fi n -i) presented by the generalized rising factorials. The 
corresponding delta operator is q{D x ) where q(t) 6 k[[t]] is given by q{t) ~ 

Proof. Evaluating q{D x ) on p n (x), we find 
e- x '» D - p n [x + -\.n)-p n {x) 

-Pn(X) 



— l.fi —1.11 

p»(-i.m) 



- X] ( n )pn-i(x)(-l.H + Pl)(-l.fl + pi +M2) •■ " (-1./U + M1 H h/Uj-l) 



*>1 



p„_i(a;) • 1 . 



Here the last line follows from the substitution lemma by replacing —l.fi + ^i 
with e = 0. □ 



Just as the transfer formula arose from the generalized Abel presentation, the 
Rodrigucs formula arises from the presentation by generalized rising factorials. 
Preserve the notation from the preceding proposition. We have 

p n (x) ~ x(x + pi)^ (x + m-\ \- Pn-l) 

~ x-e^ D *(x(x + {i 2 )---(x + p2^ hMn-2)) 

~ x ■ e IJ,Dx p n -i{x) 

- x-(q'(D)) p n _i(x), 



where the last line follows since q'(t) ~ e _1,Alt . 

We close our consideration of generalized rising factorials by observing that 

the presentation p n (x) ~ x(x + pi) ■ ■ ■ (x + fii H h/i n -i) immediately suggests 

a combinatorial interpretation along the lines of reluctant functions. 

Proposition 19 Let x be a nonnegative integer and let fi be an umbra such 
that \i % ~ rrii where each mi is a nonnegative integer. The value of p n (x) ~ 

x(x + fix) ■ ■ ■ (x + fi\ H h Hn-i) is the number of labeled forests on n vertices 

where each tree is assigned one of x colors, each vertex of outdegree j is assigned 
one of m,j colors, and where each parent vertex has a smaller label than each of 
its children. 

Proof. It suffices to observe that we can construct such a tree by specifying a 
function mapping each vertex to its parent and then choosing colors. Choos- 
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ing jij from the ith multiplicand corresponds to requiring the function to map 
vertex i to vertex j. Choosing x of course indicates that vertex i is a root. □ 



6 Presentations for Sheffer sequences. 

Recall that a sequence of polynomials Sq (x) , Si (x) , S2 (x) , • ■ ■ with s n (x) of degree 
n is said to be a Sheffer sequence with respect to a delta operator Q (or with 
respect to the associated sequence of binomial type) when Qs n (x) = ns„_i(x) 
for all n > 0. We will call a Sheffer sequence normalized when si(x) is monic. 

Let Q = f(D x ) be a delta operator. Let po (x) , p\ (x) , ... be the associated 
sequence of binomial type. Since Q and e" D:c commute, we have Qp n (x + 
0) ~ np n ^i(x + 0) thus p n (x + (3) is a Sheffer sequence for Q. Since Q/D is 
invertible, s„(0) and Qs n (x) determine s n (x). Hence si(0), S2(0), . . . determines 
any Sheffer sequence with respect to Q. But because the the Pi{x) have different 
degrees, any such sequence arises from the umbral presentation s n (x) ~ p n (x + 
0) for suitable choice of umbra /3. We have proved the following. 

Proposition 20 If po(x) , pi(x) , . . . is a sequence of binomial type, then all asso- 
ciated Sheffer sequences are umbrally presented by the sequence po{x + 0) , p\{x + 
0), P2(x + f$), ... for some a G k and some umbra D 

The following is immediate from the preceding and the presentation results 
in the preceding section. 

Corollary 21 Let sq(x), Si(x), . . . be a sequence of polynomials in k[x] and 
let po (x) , pi (x) , . . . in k[x] be a sequence of binomial type. The following are 
equivalent: 

1. The sequence So(x), Si(x), ...is Sheffer with respect to Po(x),pi(x), . . .. 

2. There exist umbrae (3,j such thatp n {x) ~ (x.7)™ and s n (x) ~ ((x+/3).7) n . 

5. There exist umbrae 0a such that p n (x) ~ x(x + n.a)"^ 1 a?i(i s n (x) ~ 
(x + /3)(x + /3 + n.a)"- 1 . 

^. There exist umbrae /i smc/i thatp n (x) ~ x(x+/ii) • • ■ (x+/ii+- • '+/x n _i) 
and s„(x) ~ (x + /3)(x + /ii) • • • (x + + fj.% H h /-t„-i). 



Since (x + t/ + /3).7 = x.7 + {y + f3.j), the standard expansion result for 
Sheffer sequences, 



□ 




follows immediately from part of this corollary. 
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7 Multiplicative sequences 



We present our final results as easy applications of the preceding constructions. 
We have relied on the umbral relation (x + y).7 = X.7 + 2/.7 to find a presentation 
for sequences of binomial type. It is natural to ask what happens if we replace x 
and y themselves by umbrae. Fix umbrae a, (3 such that a n ~ a n and (3 n ~ b n . 
We will consider sequences of polynomials in multiple variables 0,1,03,... and 
61,62, ■ ■ ■• 

Fix an umbra 7 and define K m (ai, 02, . . .) — (0.7)™. By equation |^ 
K m (ai, a<2,, . . .) has degree m, is linear in the variables a, and only depends 
on the variables ax, ... , a m . 

^From the point of view of generating functions, the fairly simple umbral 
relation (a + 0) .gamma = a. 7 + /3«7 now becomes: 

Proposition 22 Define polynomials K m £ k[ai, . . . , a TO ] 6?/ -K" m (ai, 02, ... , a m ) : 
(a. 7 ) m . // 

Z^ I z^ \ ( z^ 

fc>o y>o / y>o J 



fe>0 " \i>o J \j>o J 

We generalize this "nice" behavior of generating functions under multi- 
plication as follows. Suppose that Kq,Ki,... is a sequence of polynomials 

in the variables t\,t2, For r £ k[X], denote by K m (r) the evaluation, 

if m (E[r], E[r 2 ], E[r 3 ], . . .), of the polynomial K m . Define the sequence Kq, K\, . . . 
to be multiplicative if whenever a, /3, x are umbrae such that a + (3 = x? we 
have p + cr = r where p m ~ K m (a), a m ~ /£T m (/3), and r m ~ ^(7). 

Since n.(a + 0) = n.a + n./3, the sequence of polynomials K rn in 01, ai, • • • 
defined by K m (a\, . . . , a rn ) ~ (n.a) m is multiplicative . 

These constructions can be generalized as follows. 

Proposition 23 Let a be an umbra with a % ~ aj. Let I be a positive integer 
and 7 and umbra. For C\, . . . , cj £ k, define the polynomial K m (ai, . . . , a m ) m 
k[ai, . . . , a m ] by K m (ai, . . . , a m ) ~ £ i=1 «.(cia).7) m . The sequence 1, Ki,K 2 , ■ ■ ■ 
is multiplicative and K m has total degree m when aj is given degree i. 

7/7 = 1, then K m is homogeneous in the above grading. D 

If Kq, Ki, ... is multiplicative and homogeneous, then Lq, Li, ... where Li = 
Ki/i\ is a m-sequence or multiplicative sequences in the sense of Hirzebruch [Q, 
namely if 
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then 



y^^L k (c 1 ,c 2 ,...,c k )z k 



fc>0 



L l (a 1 ,a 2 , . . . , ai)z % ^ L j (b 1 ,b 2 , ■ ■■ , bjjz 3 
i>0 I \j>o 



8 Open problems 

We close with a brief list of open problems and areas for future work. 

1. Generalize Theorems 14 and [l?] and their corollaries by finding umbral 
presentations corresponding to other well-known sequences of binomial 
type. 

2. Determine which sequences of binomial type over the integers and which 
sequences of integral type (see [BBN]) may be presented in the above 
fashions. Find general presentation formulae for these situations. 

3. Find umbral presentation theorems for the various generalizations of the 
umbral calculus (see for instance [LR, O). 



4. Give conditions for a sequence of binomial type to be presentable as x.G 
where G is a random variable rather than an arbitrary umbra. 
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